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LS 5039 is one of the four TeV emitting X-ray binaries detected up to now. The powering source 
of its multi-wavelength emission can be accretion in a microquasar scenario or wind interaction 
in a young non-accreting pulsar scenario. These two scenarios predict different morphologic and 
peak position changes along the orbital cycle of 3.9 days, which can be tested at milliarcsecond 
scales using VLBI techniques. Here we present a campaign of 5 GHz VLBA observations con- 
ducted in June 2000 (2 runs five days apart). The results show a core component with a constant 
flux density, and a fast change in the morphology and the position angle of the elongated extended 
emission, but maintaining a stable flux density. These results are difficult to fit comfortably within 
a microquasar scenario, whereas they appear to be compatible with the predicted behavior for a 
non-accreting pulsar. 
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1. Introduction 

The very high energy emission of gamma-ray binaries (see [|l], |2|, ^ 0]) is basically interpreted 
as the result of inverse Compton upscattering of stellar UV photons by relativistic electrons. Two 
excluding scenarios have been proposed to explain the acceleration mechanism that powers the 
electrons. In the first one electrons are accelerated in the jets of a microquasar powered by accretion 
(see [|], and In the second one they are accelerated in the shock between the relativistic 



wind of a non-accreting pulsar and the wind of the stellar companion (see []10|], ||8|], and [|22|]). 

LS 5039 is a gamma-ray binary situated at 2.5 ±0.5 kpc, composed of a bright young star and 
a compact object. The mass of the compact object depends on the inclination of the orbit, which 
can be between 11 and 75°, and is in the range 1.5-8 Mq. The compact object can thus be a black 
hole or a neutron star ||6ll. 



2. Testing the possible scenarios at mas scales and previous radio observations 

The expected behavior of the radio emission at mas scales is different in each scenario. In the 
microquasar scenario we have a central core with extended jet-like radio emission that displays a 



flux and distance asymmetry produced by projection effects and the Doppler boosting [13], 



The direction of the jets can display long-term precession. The jets can be bent or disrupted due 



to the interaction with the dense stellar wind of the bright companion []18|]. In the non-accreting 
pulsar scenario the shocked material is contained in a 'bow' shaped nebula extending away from 
the stellar companion that follows an elliptical path during the orbital cycle. At mas scales, it is 
expected that the direction of the extended emission changes with the pulsar's orbital motion. The 
peak of the radio emission, located at a few AU behind the pulsar, should follow an elliptic orbit 
~10 times bigger than the size of the orbit of the system 

The unresolved radio emission of LS 5039 is persistent, non-thermal, and variable, although 
no strong radio outburst or periodic variability have been detected so far ijl^, The source spectral 



index is a = — 0.46±0.01 [11]. No radio pulses have been detected at 1.4 GHz [14]. However, 
free-free absorption with the stellar wind may prevent the detection of pulsations, and observations 
at higher frequencies are required. Observations at 5 GHz were performed using the VLBA and 



the VLA in its phased array mode in 1999 [|15p. The orbital phase of the system was in the range 
0.12-0.15. The final synthesis map shows two-sided extended emission emerging from a central 
core with a flux density of 12.8 mJy, and aligned towards a direction with Position Angle (PA) 
= 125°, measured from north to east. The source extends over 6 mas on the plane of the sky. The 
morphology of the source was studied at larger scales with the EVN and MERLIN, also showing 
bipolar extended emission emanating from a central core (see Paredes et al. [[T^). 



3. VLBA observations in June 2000 



We observed LS 5039 with the Very Long Baseline Array (VLBA) and the Very Large Array 
(VLA), of the National Radio Astronomy Observatory (NRAO), at 5 GHz frequency on 2000 June 
3 and 8. The two 8-hour observing sessions, hereafter run A and run B, took place at orbital phases 
in the range 0.71-0.79 and 0.43-0.51, respectively. The observations were performed switching 
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between the target source LS 5039 and the phase reference caUbrator J1825— 1718. The correlation 
position of J1825— 1718 was shifted by Aa = +57.0 mas and A5 = +21.4 mas from the accurate 
position provided by later observations from the joint NASA/USNO/NRAO geodetic/astrometric 
program. In the case of LS 5039, due to its proper motion, the source was found to be Aa = 
+7.6 mas and Ad = —16.2 mas away from the correlated position. See details in [|20|]. 

The post-correlation data reduction was performed using the Astronomical Image Processing 
System (AIPS). The VLA position had to be corrected by 0.43 m, according to later geodetic 
measurements. The positions of J1825— 1718 and LS 5039 were corrected using the task CLCOR. 
As recommended for phase-referencing experiments, we applied ionospheric and Earth Orientation 
Parameters corrections to the visibility data using CLCOR. The fringe fitting (FRING) of the residual 
delays and fringe rates was performed for all the sources. Fringes for 15 and 25% of the baselines 
were missing for the target source LS 5039 and for the astrometric check source J1837— 1532, 
respectively. Typical data inspection and flagging were performed. An independent reduction of 
the VLA data was performed using standard procedures within MPS. 

4. Results 

The VLA data correlated as an independent interferometer were compatible with a point-like 
source for the obtained synthesized beam of 5.5" x 3.7" in PA ~ —4°. We measured the flux 
density of the source every 30 minutes and obtained a mean of 29.4 mJy with a standard deviation 
of a=l.l mJy during the 8 hours of run A, and 28.4 mJy with a=0.7 mJy for run B. 

We show the final VLBA-i-phased VLA self-calibrated images in Fig. |. The image obtained 
for run A displays a central core and bipolar and nearly symmetric extended emission with PA ~ 
116 + 2°, with the brightest component towards the south-east. The image is similar to the one 
obtained with the same array in 1999 May 8, corresponding to orbital phases 0.12-0.15, which 
showed a slightly more asymmetric extended emission in PA~ 125° (see JTs]]). In contrast, the 
image obtained for run B displays a core and a bipolar but clearly asymmetric structure with PA 
~ 128 + 2° and the brightest component towards the north-west. The total and peak flux densities 
fitted to the data are shown in Table |l]. See also [pO|]. 

We splitted each run into 4-hour data sets, and self-calibrated each block. No significant 
morphological differences are measured between the two halves in any of the two runs. The peak 
position of the component SEl with respect to Corel is stable in 4 hours within the errors (a^ = 
0.31 mas, Og = 0.62 mas). For run B the distance between Core2 and NW2 is also stable in 4 
hours within the errors (a^ = 0.55 mas, Og = 0.25 mas). 

The measured position of LS 5039 for run A is OCj20oo.o = 1 8*^26™ 15?05653+0. 15 mas, and 
5j20oo.o = — 14°50'54''2564+1.5 mas. For checking purposes we measured the position of data 
blocks of 4 and 1-hour lengths. The peak position of LS 5039 appears to move Aa = +0.1 + 
0.1 mas, A5 = +2.8 + 0.2 mas between the two 4-hour blocks (slightly more in the 1-hour blocks). 
The observed direction has a PA of 1.7 + 1.7°, which is the same as the line joining the positions 
of the source and J1825— 1718. The expected error for differential astrometry is given by the 
separation d in degrees from the phase-reference source and the offset A of its correlated position, 
according to A x (J/l 80) x 7i (see [^]). Plugging our offset of 60.9 mas and a distance of 2?47 we 
obtain an error of 2.6 mas, very similar to the observed displacement of LS 5039. Similar results 
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Figure 1: VLBA+phased VLA self-calibrated images of LS 5039 at 5 GHz obtained on 2000 June 3 (left) 
and 8 (right). North is up and east is to the left. Axes units are in mas, and the (0,0) position corresponds 
to the source peak in each image. The convolving beam, plotted in the lower left comer, has a size of 
3.4x 1.2 mas in PA of 0°. The first contour corresponds to 5 times the rms noise of the image (0.08 and 
0.11 mJy beam^' for run A and B, respectively), while consecutive ones scale with 2'/^. The dates and 
orbital phases are quoted in the images. There is extended radio emission that appears nearly symmetric for 
run A and clearly asymmetric for run B, with a small change of ~ 12° in its position angle. 



Table 1: Parameters of the Gaussian components fitted to the data. Columns 3 and 4 list the peak and 
integrated flux densities of each component. Columns 5 to 8 list the polar and Cartesian coordinates of the 
components with respect to the peak position. The PA is positive from north to east. 



Run Comp. 


Peak 55 ghz 


S5 GHz 


r 


PA 


Aa 


^5 




[mJy beam^'] 


[mJy] 


[mas] 


n 


[mas] 


[mas] 


A Corel 


10.54 ± 0.08 


20.0 ± 0.2 










SEl 


1.11 ±0.08 


2.6 ± 0.2 


3.67 ± 0.08 


115.9 ± 1.7 


3.30 ± 0.07 


-1.60 ±0.12 


NWl 


0.88 ± 0.08 


1.5 ±0.2 


3.29 ± 0.09 


-63 ±2 


-2.92 ± 0.08 


1.52 ±0.14 


B Core2 


10.45 ±0.11 


17.6 ± 0.3 










SE2 


0.75 ±0.11 


1.8 ±0.4 


2.8 ±0.2 


129 ±5 


2.17 ±0.13 


-1.8 ±0.3 


NW2 


2.22 ±0.11 


3.9 ± 0.3 


2.94 ± 0.06 


-52.2 ± 1.4 


-2.32 ± 0.04 


1.80 ± 0.09 



are obtained when splitting the data from the astrometric check source J1837— 1532. Therefore, 
these secular motions appear to be purely instrumental. The errors assigned to the coordinates of 
LS 5039 quoted above are half of the total secular motion measured in 1-hour blocks. 

The phase-referenced image of run B reveals a double source with similar peak flux densities of 
3.6 ±0.1 mJy beam^', much lower than expected, and are separated 2.5 ±0.1 mas in PA of 87 ±5°. 
Moreover, 4- and 1-hour blocks reveal a fading of the western component and a brightening of the 
eastern one along the run, as well as a similar secular motion as in run A. Tropospheric errors, 
which affect the phase-referenced image and cannot be accounted for, can easily split Core2 into 
the observed double source. The precise position of the peak of LS 5039 cannot be measured in 
run B. The relative astrometry between runs in these observations is therefore not useful. 
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5. Discussion 



The observations of LS 5039 reported here, obtained with the VLBA on two runs separated by 
5 days, show a changing morphology at mas scales. In both runs there is a core component with 
a constant flux density within errors, and elongated emission with a PA that changes by 12 ± 3° 
between both runs. The brightest component in run A is towards south-east, and in run B towards 
north-west. The source is nearly symmetric in run A and asymmetric in run B (see Fig. |l|). The 
relative astrometry obtained with current data is not considered in this discussion. 

In the microquasar scenario, and assuming ballistic motions of adiabatically expanding plasma 
clouds without shocks [13], the morphology of run A can be interpreted as a double-sided jet 
emanating from a central core with the southeast component as the approaching one, whereas in 
run B the northwest component would be the approaching one. We can compute the projected bulk 
velocity of the jets from the distance and flux asymmetry. The results are compatible with a mildly 
relativistic jet with a projected bulk velocity of 0.1c. However, the distances from Core2 to the 
components NW2 and SE2 are very similar and do not imply any significant relativistic motion. 

The lack of proper motions (see upper limits in Sect. ^ sets an upper limit to the projected bulk 
velocities of the components. For the measured flux asymmetries, the SEl jet should be pointing 
at 6 < 48°, and the NW2 jet at < 45°. In this context, a fast jet precession is needed to explain 
this behavior, and the PA of the jet should vary considerably, in contrast to the small range covered 
by all observed values at mas scale, between 115 and 140° ||15, 16]. See discussion in |2C]. 

Alternatively, the morphology detected in run B could be the result of a discrete ejection where 
Core2 is the approaching component and NW2 the receding one, while there is no radio emission at 
the origin of the ejection. However, large X-ray and radio flux density variations are expected dur- 
ing discrete ejections (see and references therein), while the peak and total radio flux densities 
of LS 5039 are strikingly constant (see also [19] and [^) and there is no evidence of an X-ray flare 
in 11.5 years of RXTE/ ASM data. Although not yet explored in detail, the morphology changes 
can be due to the interaction between the jets and a clumpy or dense stellar wind []T^]. 

In the non-accreting pulsar scenario, the different morphologies we have detected at different 
orbital phases could be due to the change of the relative positions between the pulsar and the com- 
panion star along the orbit (see details in [^]). Observations at different orbital phases have always 
revealed a very similar PA for the extended emission, which would correspond to an inclination of 
the orbit as seen nearly edge on (90°). On the other hand, the absence of X-ray eclipses places an 
upper limit of / < 75° in this scenario. Therefore, these two restrictions imply an inclination angle 
that should be close to the upper limit of 75°. 

In conclusion, a simple and shockless microquasar scenario cannot easily explain the observed 
changes in morphology. On the other hand, an interpretation within the young non-accreting pulsar 
scenario requires the inclination of the binary system to be very close to the upper limit imposed 
by the absence of X-ray eclipses. Precise phase-referenced VLBI observations covering a whole 
orbital cycle are necessary to trace possible periodic displacements of the peak position, expected 
in this last scenario, and to obtain morphological information along the orbit. These will ultimately 
reveal the nature of the powering source in this gamma-ray binary. 
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